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Uniform Asymptotic Expansions for Weber 
Parabolic Cylinder Functions of Large Orders 

F. W. J. Olver* 

(April 6, 1959) 
Expansions of solutions of the differential equation 

d 2 W 



dt 2=^\t 2 -l)W, 

are sought for large values of |/*|, which are uniformly valid with respect to arg n and un- 
restricted values of the complex variable t. Two types of expansion are found. Those of 
the first type are in terms of elementary functions and are valid outside the neighborhoods 
of the points t= ± 1, the turning points of the differential equation. The second are in 
terms of Airy functions and hold in unbounded regions containing one of the turning points. 
The special forms of the expansions when the variables are real are considered in detail, 
and asymptotic expansions for the zeros of solutions of the differential equation are found 
by reversion. Numerical examples are included. 
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Part 1. Introduction 

1. Introduction and General Summary 

The purpose of this paper is to determine asymptotic expansions of solutions of the differ- 
ential equation 



5=^ 2 -^ 



(i.i) 



for large values of the complex parameter jjl which are uniformly valid with respect to arg /x 
and the complex variable /. 

The immediate reason for this investigation is that the solutions of eq (1.1) are basic 
functions in the theory of the asymptotic solution of the differential equation 



^={Mt)+i(t)}w 



■ i.-j 



for large values of [/x] in a region in which p(t) has two simple zeros, just as Airy functions are 
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basic functions in the theory of eq (1.2) in a region containing one simple zero of p(t). A 
knowledge of the uniform asymptotic behavior of the solutions of (1.1) is required in developing 
the asymptotic theory of (1.2). 

Equation (1.1) is of interest, however, in other connections, for example in the study of 
wave motion inside or outside parabolic cylinders, and wave propagation in an inhomogeneous 
atmosphere. Hermite polynomials of order n are expressible in terms of the solutions of the 
equation when iJL 2 =2n J r 1. Accordingly, for the sake of completeness many results are derived 
in this paper which go beyond those needed for the immediate purpose stated in the preceding 
paragraph. 

The desired asymptotic expansions are obtained by application of the general theory of 
the asymptotic solution of linear differential equations of the second order as developed by the 
present writer in [9]. 1 From the standpoint of the general theory, eq (1.1) is characterized 
by having turning points at t=±l. Asymptotic expansions in terms of elementary functions 
can be found in any region bounded away from these points and expansions of this kind are 
obtained in part 2 by application of theorem A of [9]. 

In order to determine the behavior near the turning points, expansions in terms of Airy 
functions need to be used. In part 3 these expansions are obtained by application of theorem 
B of [9] ; they are Valid in unbounded ^-regions which include one (but not both) of the turning 
points. All combinations of fi and t are covered. 

When the variables are real, eq (1.1) has important special solutions. Expansions of 
these solutions are derived in part 4 from the results of parts 2 and 3. Also included in part 4 
are numerical examples illustrating the powerful nature of the expansions. 

In part 5 uniform asymptotic expansions for the zeros and associated values are obtained 
by reversion and illustrated by numerical examples. Much of the analysis of this part and 
some in earlier parts is similar to that used in [8] in the theory of Bessel functions of large 
order, and advantage is taken of this analogy whenever possible. 

Investigations of the asymptotic solution of eq (1.1) for large \fx\ have been made by 
several writers. Their work on expansions in terms of elementary functions is described in 
section 6 of part 2, and on expansions in terms of Airy functions in section 10 of part 3. At 
this stage it suffices to remark that the expansions obtained in this paper include practically 
all of the earlier results and extend them in various ways. 

2. Relevant Properties of Parabolic Cylinder Functions and the Gamma Function 

A comprehensive account of the properties of parabolic cylinder functions is given by 
Miller in the Introduction to [4]. For the purpose of reference we collect here the properties 
which will be required in this paper. We use Miller's notation rather than that of Whittaker 
[13] and adopt Miller's choice of solutions when the variables are real. 

The standard form of differential equation for the parabolic cylinder functions is 

The principal solution U(a,z) is determined by the condition 

U(a,z)~z~ a -h-i z2 as z->+ ». ( 2 . 2 ) 

It is an integral (entire) function of z and an integral function of a. In terms of Whittaker's 
notation D n (z) for parabolic cylinder functions and W k>m (z) for the confluent hyp ergeome trie 
function, we have 

U(a,z)=D-a-i(z)=2-i a z~*W-ia,-i Q z 2 \ (2.3) 



i Figures in brackets indicate the literature references at the end of this paper. 
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Other solutions of (2.1) are U(a, — z), U(—a,iz) and U(—a,—iz). The connection 
formulas are 

U(-a,±iz) = (2w)-*T (^+a){e- iir ^ a ~^ U(a y ±z)+e i *& a ~V U(a,Tz)}, (2.4) 

6 7 (a,±£) = (27r)--*rQ-aW (2.5) 

For fixed a and large \z\ 

U(a,z) = z- a -h-* z2 {l + 0(\z\- 2 )} (|args|< |ir— € ), (2.6) 

where 2 the is uniform with respect to arg 0. Here and elsewhere e denotes an arbitrary 
positive number which is independent of all other variables. From (2.4) and (2.6), we obtain 

^(a- S )=e- i ' to H) 2 -«-V^ 2 {l + 0(| 2 r 2 )}+-^-3''-i e ^{l + ()(| 2 r 2 )}, (2.7) 

valid when — |7r+e<arg z<\tt— e. The same result with the sign of i changed holds when 
— i7r + e<arg z<%ir—e; the apparent discrepancy between the formulas in their common 
region of validity is merely an example of the Stokes phenomenon. 

Let the 2-plane be divided into four sectors M () , M b M 2 , and M 3 , defined by |arg z\ <|7r, 
i7r<arg z<\-k, I arg {—z)\ <\ty and — |7r<arg z< — \ ir, respectively. Then eq (2.(5) shows 
that for large \z\, [J(a,z) is exponentially small in M and exponentially large in Mi and M 3 . 
Equation (2.7) shows that U(a,z) is also exponentially large in M 2 , unless a has one of the 
values — I, — f, — |, . . . , 111 which case it is exponentially small. In fact in this event 

U{-n-\-z) = (-YU{-n-^z) (n=0,l,2, . . . )• (2.8) 

The solutions U(a,±z), U(—a,±iz) clearly form a numerically satisfactory set in the 
complex plane, because in each sector M ; one of the solutions is exponentially small and at 
least two others exponentially large. Accordingly, we may expect that a knowledge of the 
asymptotic behavior of U(a,z) for large \a\ in the sector |arg z\<%ir } or even the half-plane 
|arg z\<^w, will enable 4 us to determine the behavior of any solution of (2.1) over the whole 
£-plane by use of connection formulas. 

Solutions when the variables are real. There are two forms of eq (2.1) of importance when 
the variables are real, namely 



and 



d 2 w 
dx 2 



=(a—-x 2 Jw, (2.10) 



in which a and x are real. Equation (2.10) is obtained from (2.1) by setting z=xe~\ iri and 
replacing a by ia. 

Keal solutions of (2.9) are U(a,x) and U(a,—x) but they do not form a satisfactory pair 
for all values of a; in particular they are linearly dependent when a=— ^, — f, — f, .... 
The fundamental solutions are taken to be U(a,x) and V(a,x), where 



V(a,x) = (1/tt) t( -+a Vsin ira- U(a,x) + U(a, 



-x)}. (2.11) 



2 Miller [4, p. 27], who is more concerned with asymptotic expansions fulfilling the uniqueness conditions of Watson [llj than the Poincare 
condition, gives the range |arg z|<§tt. The validity, in PoincarS's sense, for the wider range is established in [14, p. 347]. 
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For all values of a, U(a,x) and V(a,x) are linearly independent. When a is negative it is often 
convenient to use in place of V(a,x) its multiple U(a,x), given by 

U(a,x) = Ty y— a )V(a,x) = t&n 7ra-C7(a,x)+sec Tra-U(a,—x). (2.12) 

When a+i is a positive integer U(a,x) is infinite. 
Wronskian relations for these solutions are 

W{U{a,x)y{a,x) } = (2/t)*, i^{U(a,x),U(a,x) } = (2/x)*r (5-A (2-13) 

The Hermite polynomials 

ff w (z) = (-)V 2 ^<T* 2 (2.14) 

are related to the function U(a,x) by 

H m (a!)=2*»«* a 'i7(-n-|,ah/2). (2.15) 

For eq (2.10), the fundamental pair of real solutions is taken to be W(a,x) and W(a, — x), 
denned by 

W(a,x) = (2k)^e^ a Ue{e i{ ^^ ir) U(ia ) xe-i 7ri )} y (2.16) 

W(a-x) = (2/k)hi 7ra Im{e i{ ^ + ^ ) U(ia,xe-^ i )}. (2.17) 

where 

fc=(l+e 2rt )*— e", l/^(l-f^ fl )Hr <fe=arg r Q+m\ (2.18) 

The value of arg rQ+ia) here is not the principal one but is determined by the conditions that 
it is a continuous function of a and equals a In a— a J rO(a~ 1 ) as a — >+ 00 . 
The Wronskian relation is 

7T{W r (a,a;) l W(a,-a;)} = l. (2.19) 

Gamma Junction expansions. We shall need the asymptotic expansion of the function 
t(--\-z ) and its reciprocal for large \z\. From Stirling's series [5, chapter 9], we have 

In rQ+A^g In 2t+z In z-z+Z(z) } (2.20) 

where 

Z(z)~± J*® ^4l+oLklL js+ ■ . • (|arg S |<x- e ) (2.21) 



uniformly with respect to arg 2, B 2s denoting the Bernoulli polynomial of degree 2s. Hence 

r(l+z)~(27rye-*z*j:^ (largsl^x-e), (2.22) 

where the constants are determined by the asymptotic identity 

S5~exp{Z( 2 )}, (2.23) 

or equivalently, 

§>-{SW}{£S}- (2.24) 
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Numerical calculation yields 



_i l l l,)(,;l _ 4027 , OOR . 

To-1, ^1--^ ^-Ii52 ? 73 -4T4720' 74 ~ ~398 13120* (2 ' 25J 

Since Z(z) is an odd function of z, it follows from (2.23) that 

{SifMs*-)'?}- 1 ' (2 - 26) 

identically, and hence that 

-^— -|^r 2 ± (-)' g (|arg .| <»-«). (2.27) 

r(J+s) (2?r)2 5 =o ^ 

Part 2. Expansions in Terms of Elementary Functions 

3. Asymptotic Solutions of the Differential Equation 

In order to bring the differential eq (1.1) into the standard form required for the applica- 
tion of theorem A of [9], we take new variables £=£(£) and W (not t he \\'(a,x) of sec. 2), defined 

by 

(§) 2 =< 2 -i, w=(f-iyw (3.i) 

(cf. [7] sec. 2). Equation (1.1) then becomes 

dW 
where 



[S+^mW, (3.2) 



«r^) = (« 2 -l)- i | S («»-l)*< 4^=^5- (3-3) 

For convenience we prescribe that £=0 when £=1. From the first of eqs (3.1) we then 
obtain 

6= f ' (* s — 1)* il*=* *(**— 1)*— * ln{ *+ (* 2 — 1)* } , (3.4) 

and on expansion 

6-if-ilnM-i+i g 1 ^4,i ) -,gy 1) (l«l>D. (3-5) 

Temporarily we suppose that the branches of the many-valued functions occurring in (3.4) and 
(3.5) are positive when £>1 and determined by continuity elsewhere. The relation 

*«**•)=*(*) =F*tr (3.6) 

follows immediately from (3.5) when |£|>1, and thence by analytic continuation throughout 
the £-plane cut along the join of t= — l and t=-\-l. 

The mapping of the Z-plane, cut along the real axis from — a> to + 1 , is illustrated in figure 1 . 
The center diagram corresponds to the upper half of the 2-plane and the right-band diagram 
to the lower half. They are of course conjugate to each other. Features of the mapping are 
that the upper side of the real Z-axis is mapped on three straight lines BA (Im £=0, Re £>0), 
AE (Re £=0, 0>lm £>— ^7r), ED (Im %=— iw, Re £>0), and that the positive imaginary 
*-axis is mapped on the straight line OC (Im £— — ^7r, Re £<0). 

Next, we make the substitutions (cf. [9], sec. 1) 

u= | M | 2 , 0=arg /*, Z=e- 2id x (3.7) 
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£-plane (i) £-plane (ii) 

Figure 1. Mapping of the function £(t). 

(the x here is not to be confused with the real variable of sec. 2). Then (3.2) becomes 

^={u 2 +e- Ai9 &'(e- 2iB x)}W 9 (3.J 

which has the form of eq (2.1) of [9]. The relationship between a and t is given by 



x=x(t)=e 2ie £(t)=e 2i$ 



(t 2 -l)klt; 



(3.9) 



the £-map of the £-plane is obtained from the £-map by rotation through an angle 20. The 
function jFis expressed in terms of t by (3.3). For large \x\ we deduce from (3.3) and (3.5) that 



4 16a: 2 



(3.10) 



uniformly with respect to bounded arg x and bounded 0. 

In consequence of the last result we could, if we wished, take the ar-domain D(0) of [9], 
section 2 to be the map of the £-plane cut along the real axis from — oo to +1 with the circles 
|£±1|<<5' removed, 5' being a fixed small positive number. The preliminary condition (2.2) 
of [9] would then be satisfied with «r=l. Greater regions of validity in the desired expansions 
will result however, if we arrange that the boundaries of D(0) are parallel to the positive imag- 
inary a;-axis where possible. Accordingly, we shall now investigate what we shall call the 
principal curves in the £-plane, namely the £-maps of the straight lines in the £-plane passing 
through A (£=0) and E (Z=Tiiir) and inclined at an angle —20 to the imaginary £-axis. 
In other words, the principal curves are the level curves of the function e xit) which pass through 
t=±l. 

Near t—l we find from (3.9) that 



a?(0^2V*(«-l)*. 



(3.11) 



Hence three principal curves emanate from t=l, and have as tangents the rays 



arg(*-l)=±ir-|tf l *■- 1*. 


i 


Equation (3.6) shows that three principal curves also emanate from t= — l. 
For large \t\ we have from (3.5) 




*(*)~!«*V. 


(3.12) 
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Hence the six principal curves are asymptotic to the four rays arg t=±\it — $, ±f7r— 9. 

Figure 2 illustrates the principal curves when < 0<>. Those through A are the continuous 
curves AP, AQ, AS when O<0< %t, and AP, AQ, AR when i7r<0<7r. When 0=0, AQ 
degenerates into the curve EQ plus the join AE, and when 0=|-7r, AS, AR degenerate into 
AE+ES, AE+ER, respectively. The principal curves through E are the images in the origin 
of those through A; this can be verified by using the continuation formula (3.6). 

We can now define the branch of the many-valued function £(t) which will be the most 
convenient for our purposes. When 0<#<V we cut the /-plane along the principal curves 
AQ and ER and select the branch of £(£) which is continuous in the cut plane and for which 



1. 



1. 



£(0) = >- (O<0<tt), {(0-0i)=-rwr (0=0) 



(3.13) 




(a)0=O 




[b)0<6<±Tr 




(c)e--±Tr 




(d)j-7r<e<ir 



519835—59- 



Figure 2. t-plane: principal curves and domains S(0). 
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When — 7r<0<O the cuts are taken to be the conjugates of those for — 0, and the branch is 
fixed by 

«0)=-|6r (-tt<0<O), «(0+0i)=:-jiT (0=0). (3.14) 

It is evident from (3.6) that (3.13) and (3.14) are consistent when 0=0. Henceforth we under- 
stand £ and £(t) to denote this branch. For a given value of 6 (— 7r<0<7r), £(f) is a single- 
valued function of t defined everywhere in the £-plane except when 0=0 and t lies in the shaded 
region to the left of ER + EQ in figure 2(a). When £>1, £(t) is positive provided that |0|<Cf?r; 
this result is false if %Tr<^\d\<^ir for the cut AQ passes across the real axis as passes through fx. 

Let x=x(t) now denote the branch e 2ie %(t). We then take the ^-domain D(0) of [9], 
section 2 to be the map of the function x(t) from which circles surrounding the points A and E 
and of radius ^-5 have been removed, <5 being a fixed small positive number. 

The map of x(t) has two sheets when 0=0 and three sheets when 0< |0| <V- It is illustrated 
in figure 3 for 0=0 and O<0<-^7r. The passage from one sheet to another takes place across 
the positive real axis and, when O<0<^7r, the join AE. The heavy lines correspond to the 
cuts in the £-plane and are boundaries of D(0). 




O<0< y7T(j) 



0<6< yir(ii) 

Figure 3. x-plane: domains D(0) and G(0). 
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o< e < — Trdii) 



The map of x(t) is clearly a single unbounded domain. In order to ensure that D(0) is a 
single domain we have to impose the restriction |0| <t— e, omit sheet (hi) when |0|<e and make 
8 = 8(e) small enough. Otherwise, for example, the circle center E radius %8 would intersect 
the boundary AQ. 

Next, we define the ^-domain G(0) required by [9], section 2 to be that part of the map of 
x(t) which remains after removing all points at a distance less than 8 from its boundaries. 
Sheet (hi) is again to be omitted when |0|<V The boundaries of G(0) are the dotted lines and 
circular arcs of figure 3; their £-maps are the dotted curves of figure 2. Clearly when |0| <t— e, 
G(0) is a single domain if 8= 8(e) is small enough. We readily show that conditions (i) and 
(ii) of [9], section 2 are fulfilled (cf. [9], sec. 3(i)). 

Following [9], section 2 further, we take a 2 (0) to be the point at infinity on the positive 
real axis. We then find that H 2 (0)=G(0). From theorem A it now follows that there exists a 
solution W 2 of eq (3.8) such that for large positive values of u 

W 2 ~e-™1t(-Y ^* {e s lX) (x*G(0), |*I<t-«), (3.15) 

uniformly with respect to x and 6. The coefficients are given by the recurrence relation 

j*Ui(*f*)= -\ f x J*** ( e > *)+| J e~ m ^(e- 2i9 *)^f (d, z)dx, (3.16) 

the path of integration lying in D(0). 

If we now restore the variables /* and £ by means of eqs (3.7), and let 

s/ t (() = (-y^'s/iiO, x), (3.17) 

the relations (3.15) and (3.16) become 

^~«" A E# «eS(arg M ), |ar gj u|<x-«) (3.18) 

s=() /x 

as |/x| — > oo 7 where 

V, ,,(£)=4 ^(&-\^m)s/ s (Qdi. (3.19) 

Here S(arg jjl) =S(0) denotes the £-map of G(0) ; its boundaries are the dotted curves of figure 2. 
It is of some interest to determine the asymptotic behavior of the boundaries of S(0). 
Consider, for example, the doited curve AQ^ Its equation may be written 

Z=(y-i8)e ( ^- m , (3.20) 

where the real parameter y ranges from to +°°. Reversion of (3.5) yields 

^=2{+iln8J+g+0(r 1 ln Q as |*|->». (3.21) 

When y is large, we may substitute (3.20) in (3.21) and thence obtain 



Kg ,_| ,-,+^9 !ni i , + cos28+^-f,),in2»-« + 



>my 



The corresponding equation of the principal curve AQ is of course obtained on setting 5—0 
in this result. Clearly the principal curve AQ and the dotted curve AQ! ultimately always lie 
on the same side of their common asymptote, except when 0=f tt, in which event AQ is its own 
asymptote. 
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4. Identification of Solutions 
The substitutions 

a=-\ix\ z=/**V2, (4.1) 

transform eq (2.1) into (1.1). This result and the second of eqs (3.1) show that 

TT=(^-1)*Z7(~mV«V2) (4.2) 

is a solution of (3.2). Equation (2.2) shows that when n is fixed and t->coe~ iarg ** this solution 
becomes exponentially small. The same is true of the solution W 2 introduced in section 3, 
provided that |/x| is sufficiently large and |arg fx\ <ir—e; this can be seen from (3.18) and (3.5). 
Since all solutions which are linearly independent of (4.2) are exponentially large in the same 
circumstances, we see immediately that the solution (4.2) is a multiple of W 2 . 
Hence we have 

^^MV^J-ffW^F^S^^ (*eS(arg M ), |arg M |<7r-e) (4.3) 

as |/x|— >oo , uniformly with respect to t and arg /x. The branch of (t 2 — 1)* here and elsewhere is 
understood to be continuous in S(arg /x) and to be asymptotic to the principal value of -yjt as 

t-^> oo e -i™zn t Similarly for the branch of (t 2 — 1) * used later. The function g{yi) is independent of t. 
Its asymptotic expansion for large |/z| can be found from the limiting form of the ratio of the 
two sides of the relation (4.3) as t-> oo e ~ i arg *. Thus with the aid of (2.2) and (3.5) we find that 

1 ^2*^^/r**H S % (|arg M | <*-e) (4.4) 

as |jLt|— >oo ? uniformly with respect to arg jit. Here 

g 8 = Urn ja£(£) as £-» oo <r 2 * ar g m. (4.5) 

That this limit exists and is independent of arg /x can be seen from the lemma of [9], section 5. 
The g s are not determined however until the arbitrary constants in the recurrence relation 
(3.19) have been fixed, and we now consider how to do this in the most advantageous manner. 
From (3.19), (3.3) and the first of (3.1), we obtain 

. 1 1 ds/ s , 1 f 3£ 2 +2 . . 

Taking j^!=l, we find 

* 3 -6* 



24(£ 2 -lp 



constant. (4.7) 



A convenient value for this constant is zero, then s/ x has the form of a polynomial in t divided 

by P-i)*. 

The eq (4.7) suggests the substitution 

■stf^.Cfl/P-l)*' (4.8) 

in (4.6). This leads to the recurrence relation 

(fi- l)u' 8 (t) -Zstu s (t) =r f _! («), (4.9) 

where 

8r s (0-(3/ 2 + 2)^W-12(s+l)^_ 1 (0+4(^-l)K_ 1 (0. (4.10) 
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We can satisfy this relation by taking u s (i) to be a polynomial in t of degree 3s which is an odd 
or even function of t according as s is odd or even. All the coefficients are determined auto- 
matically if s is odd, but there is a degree of freedom if s is even, since the left of eq (4.9) 
vanishes when u s is a multiple of (t 2 — l)* s . We remove this by making the coefficient of the 
highest power of t vanish. 

Thus we prescribe that u s (t) is a polynomial in t of degree 3* (s odd), 3s— 2 (s even, .s>2). 
With u (t) = l this condition determines the coefficients completely. Numerical calculation 
based on (4.9) and (4.10) yields 



u (t) = l, u l (t) = (t 3 -6t)/24, u 2 (t) = (-9t 4 +2m 2 +U5)Jll52, 
u 3 (t) = (-±0m»+18189t 7 -28287t 5 -l 51995£ 3 -2 59290£)/4 14720. 



(4.11) 



Hence 



1 9 9 1 

g 2s =0 (s=l,2,...), 00=1, 0i=2? g3= "" 2 07360 ' (4 * 12) 



An independent method of calculating the asymptotic expansion of g(ii) for large /x will be 
given in section 5. 

Expansion for the derivative. Theorem A of [9] shows that the asymptotic series (4.3) 
may be differentiated term by term with respect to t. With the aid of the first of eqs (3.1) 
we deduce that 

U'(-\»\y.m^~-^g{,x){t>-\)h-^^'^ (* e S(ar gA 0, |arg j«| <*-«) (4.13) 

as |ax| — > oo ^ uniformly with respect to t and arg fx. Here, in the notation of Miller, 

U'(a,z)=c)U(a,z)/dz. 
The coefficients are given by 

*tt)=^(*)/« a -l)*', (4.14) 

where 

v s (t)=u s (t)+^ tu s . x {t)-r s . 2 {t). (4.15) 

The first four v s {t) are 

Vo (t) = 1, vtf) = (f + 6t)/24t, v 2 (t) = (15* 4 -327* 2 - 143)/1 152, 1 



tf 3 (*) = (-4042£ 9 +18189* 7 -36387* 5 +2 38425£ 3 + 2 592900/4 14720. 
5. Use of Connection Formulas 



(4.16) 



In this section we seek expansions which hold in the complement of the domain S(arg /x), 
Since U(a,z) is an integral function of a it will clearly suffice to obtain a set of expansions which, 
together with (4.3), cover the whole £-plane (except of course the neighborhoods of the points 
£=±1) when |arg/x|<j7r. 

For convenience we introduce the notation U,- = U^(0) (j=0, 1, 2, 3, 4) for the closed 
domains into which the 2-plane is divided by the principal curves defined in section 3. Figure 4 
shows the enumeration of these regions when O<0<V. When —t<^6<0 we define 

U,(0) = U?(-0) = 0,2,4), V 1 (d) = VU-e), U 3 (0) = Uf(-0), (5.1) 

where the star denotes the conjugate region. Then the boundaries of each region \Jj(0) vary 
continuously with 6 in the interval (— t,t). 
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U 2 E\ (A U, 




(a) «=o 



(b) 0<8<-£ir 



(d)f OO 



Figure 4. t-plane: domains U;(0). 
The following relations, suggested by figure 4, are easily verified: 

Uo(0+^r)=U 3 (0), U 1 (»+|x) = UoW, U 1 (*+ir) = U 1 (fl) > 

u 3 (e+~x)=u 2 (e), u 4 (ff+gT)=u 4 w > 



(5.2) 



Uo(0+x) = U 2 (0), U I (tf+x) = U,(fl), U 2 (9+t) = UoW, 

Ug(e+») = u 1 w, u 4 (e+x)=u 4 (e). 



(5.3) 



,"A 



Figure 5. t-plane 

(o<*<^) 



C(V UV 2 ) 



The region formed, by adding to U^(0) points whose x-maps (cf. fig. 3) are at a distance 
less than 5 from the #-map of the boundary of U,(0) we shall denote by V ; -(0) or V ; . Figure 5 
illustrates for example V o (0), V 2 (0), and C{ V o (0)UV 2 (0) }, the complement of the union of V o (0) 
and V 2 (0), when O<^0<Cj7r. The boundaries are of the same character as the dotted curves 
of figure 2. 

The first connection formula we use is 



obtained from (2.5). Replacing /x by =F^m in (4.3), we have 



(5.4) 



(5.5) 
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uQfM^ygd,) -** § (-)«^ (* * s(«+±,) f 



« 2 -l)i 
and from (4.4) and (4.12), 



-2 «•+«<' 



'4*-') 



\ 3 = M / 



' 1 3 

3 1 



(5.6) 



(5.7) 



When \B\<\tt— e, we can substitute (5.5) and (5.6) in (5.4). Care is required, however, because 
at any given point in the £-plane the branches of £=£(£) and (t 2 — 1) ¥ may not be the same in 
the two formulas. 

When O<0< \-k— €, S(0— §7r) is the complement of the region enclosed by "dotted curves" 
surrounding the principal curves AS and ER of figures 2(a) and 2(b); S(0+I?O is the comple- 
ment of the region enclosed by dotted curves surrounding AP and EQ. To the right of 
AP+AS, i.e., in U o (0), the branches of £(£) and (t 2 — 1)% are the same in (5.5) and (5.6). Sub- 
stitution of these formulas in (5.4) leads to a nugatory result, however, as we would expect, 
since e^ is exponentially large here whereas U( — \ix 2 , fxt^2) is exponentially small. 

In the central region between AP + AS and EQ+ER, i.e., in C(U UU 2 ), the two branches 
of £(£) have opposite signs. Let 



£=—ir) in (5.5), £=ii? in (5.6). 
Then 77 is real and positive when — 1<7<1; from (3.4) we find 

ti= f\l-t 2 )idt=^cos-H--^t(l-t 2 )i. 

The appropriate branches of (t 2 — 1)4 and j4(£) in (5.5) are given by 

(t 2 -l)l=ei iv (l-t 2 )±, ^ s (i)=e-^u s (t)l(l-t 2 )^ 



(5.8) 



(5.9) 



(5.10) 

(cf. (4.8)), in which (1— 1 2 )± and (1 — t 2 )^ s have their principal values when — 1<*<1. For 
(5.6) we change the sign of i. 

Carrying out the substitution, we obtain 

U{— 2M , pty2) ~— — —7 cos (m ^ — **) 2-i 7T 

(1 — t 2 )* L s=o U- 



'^M%-sm(Sv-l*)± 



-^"V 



(l-i 2 ) 3s +4 



■lit) 1 

V S+2 J 



(5.11) 



as |/i|— *°°, uniformly with respect to i e C(V UV 2 ) and 6 in the interval ( — $«■+«, §«—«). Here 

?(M)~(2x)-ira+V)2^-MV^- i /i-±%^} 



,2-l" l -i«-l^^-*|± A} {l-± 5^} (I*I<*t-«), 

" Um ; J v. s=o m J -i 



V.s = 



(5.12) 
on substituting for r(J + ^ 2 ) by means of (2.22). 

The expansion (5.11) is one of the required formulas. Before leaving it we shall show that 
g(ii)=g(n) ) and in the process obtain a new derivation of the coefficients in the asymptotic 
expansion (4.4) for l/g(n). 

When e<e<\ir— e the principal curve AQ is included in the region of validity of (5.11); 
this is easily proved by considering the #-map. On this curve y 2 ^ is positive. Let us set 



„=„„=( m +i)-? 



(5.13) 
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where m is a positive integer. Then (5.11) reduces to 

Ji ( l 2 * /9 \ ( n w 2g(/x) ^ (-) s u 2 s (t m ) ( , 

where £ m is the value of £ corresponding to r]=rj m . 

Now keep /x fixed and let ra->oo. Then t m -^<*> e 1 ^*'^ , and so from (2.7) 

U (-i m 2 , ^W2)~e-^(M^V2)^^+ (2 ^ ( ^ §) «**a(rt.V2)-*^*. (5-15) 

From (3.5) and the first of (5.8) 

\ tl-\ In 2t m =-i Vm +\+0(\t m \- 2 ). (5.16) 

Hence 

e-^ H Hnt m V2)i" 2 ~i~ (_)-e^-i" 2 2-J" 2 "i m**'"*^*, 

eJ"^(M^V2)-i" 2 -'~(-) m e- i ' i+J '' J 2i'' 2 -^~ i '' 2 -^; i > 
and substituting these results in (5.15), we obtain 

( _ )mei ^ 2 -v-^-* r 1+ wy^ L- t| (M7) 

L ■*■ \2 2M ) J 

as the limiting form of the left-hand side of (5.14). 
On the right of (5.14) we have 

(l-t 2 m )i~e-^tL u 2s (t m )l(l-a)* s ->Q («>1) 
(cf. (5.10), (4.5) and (4.12)). Comparing the two sides we see that 

5M ^-i-» t -i-V-'{| + gj >^7_y' } «<«<|*-0, (6.18, 

where we use the symbol ~ to mean that the two sides of this relation have the same asymp- 
totic expansion for large |/i|. 

We now substitute (2.27), with z=— ^ jjl 2 , in (5.18) and obtain 

Although this result has been proved on the assumption that e<0<i-7r— e, comparison with 
(5.12) shows immediately that it holds in the range |0| <-| w— e, and that 

identically. Thence using (2.26) and changing the sign of n 2 , we find 

Multiplying this identity with (5.19), and comparing the result with (4.4) (remembering that 
<?2s=0 if s>l), we see that we can take 

000=000 (5-22) 
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in the expansion (5.11), as stated earlier. 

We can express the coefficients g 2s +i in terms of y s , by equating coefficients in (5.21); thus 



0i=— 7i= 



1 
: 24' 



-2^72—473=- 



2021 



2 07360 



on using (2.25). These values agree with those given in (4.12). 

The expansion (5.11) is the first of our expansions for U(— \\£, ixt^2) in regions comple- 
mentary to S(0). The corresponding expansion for the derivative is 



U'(-y\ ^2)~^2g(»)(l-t 2 )i [sin (/x 2 ^-^) S 

L *=o 



(-) S V2s(t) 

(i-* 2 )V s 



+ COS( M 2 77— \tt) 2 



(— ) s V 2s +i 



f=5 (1— t*) 8t+ * 



M 4S+2 J 



(5.23) 



valid under the same conditions as (5.11); the coefficients v s (t) are given oy (4.15) and (4.16). 

A further expansion can be found by substituting (5.5) and (5.6) in (5.4) when t lies in the 
region (U 2 ) to the left of EQ+ER. We shall not carry out this substitution however, because 
the expansion so obtained is contained in a more general result (5.28), which we shall now derive. 
As a preliminary we record the following relations in which the symbol ca has the same meaning 
as in (5.18). 



gW) : 



MW-i) 



g(v) 



(— 7r+€<argM< — c), 



(2w)ie~ iT ^ fi2 ^ _i- 1 1 

— jTTi — r^ — g(w ***) ~ 2 cos -V^(m) (|arg m| <^~ «), 



We use the connection formula 



:e ^Vm) («<arg /x<7r— e). 



(5.24) 
(5.25) 



^(^|^-M*V2)-^ ( "*^* ) ^(-|^iM^)+ (2 %l-Q) " uQf-ht&y (5.26) 

obtained from (2.4). We may substitute the expansions (4.3) and (5.5) when 

— |tt+€ <0<tt— e. 
The branches of £(t) are the same in both expansions when 

t e C(y 2 UY 3 \JY 4 ) (o<*<|A 

Z€<7(V 2 UV 3 ) (-^7T+€<0<O, \*K<e<TT-^ 

Using (5.25) we obtain immediately 

U (-i m 2 ,- W2) ~ -g-j [(sin I ,„»+* cos| v) •" 



(5.27) 



+2cosixM 2 -e" £ 2:(-) 



*=o M 



v .^(ei 



2S 



*e<7(V 2 UV 8 ) (-|ir+«<0<o)> 

<eC(V 2 UV 3 UV 4 ) (o<e<ix-A 



(5.28) 



519835—59 6 
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and 



z7 (~4 m2 ~ m<V2 ) 



^, e -¥^ 



g(p) 



(t 2 -i) 



i)«L 



■" 2| T, ^f-+e" 2( ±(-y 

s=o M s=0 



Ml) 11 
m 2s J 



< e <7(V 2 U V 3 U V 4 ) (e<6< \ A 
«6C(V 2 UV 3 ) (J x<0 <*-«)• 



(5.29) 



The last result may be expressed in a form more pertinent to the range |0|<§x by replacing 
ix by /j.e iw and using (5.24). This gives 



°U«**h&zl 



— „ 2 

s=0 M 



s=0 



j*C(«)ll 



teCiWoUYi) (-l*<e<-*y 



(5.30) 



Three further formulas, (5.28*), (5.29*), and (5.30*), can be derived immediately from 
(5.28), (5.29), and (5.30), respectively, merely by changing the sign of i. They are valid with 
the conditions 



t t OTUVjUV*) 
t e^VjUVj) 

< 6 c7(V UV 3 ) 



(o<0<|ir-A 

(-x+e<0<-;!j7r), 
(-|r<*<- € ), 

(e<0<^) 



► for (5.28*), 



for (5.29*), 



for (5.30*). 



6. Summary of Expansions in Terms of Elementary Functions 

The principal expansion of this kind is (4.3). It is valid for large \n\ uniformly with respect 
to arg n = 9 and t when — ir-\-e <0< r—e and t lies in the unbounded domain S(0), the bound- 
aries of which are the dotted curves illustrated in figure 2. The function £(t) is given by (3.4), 
that branch being chosen which is continuous in S(0) and takes the value (3.13) or (3.14) at 
the origin. The branch of (t 2 —!)* in (4.3) is continuous in S(0) and is asymptotic to the prin- 
cipal value of y/t as t-^^e~ ie Both £(<) and (t 2 — 1)* are positive when £>1 and |0|<fx. The 
coefficients «£#(£) are given by (4.8) and (4.11), the branch of {t 2 — l)^ s being defined similarly. 
The function g(n) is calculable from either of the asymptotic expansions 



—2 e „ |l + g-^|, 

—i, ,2_i — I „2 i,,2_i C 1 °° y 1 



(6.1) 
(6.2) 



both of which are uniformly valid with respect to arg /x in the interval (— ir-\-e, T—e). Here 
ft, Qz are given by (4.12); 71, Ya, 7s, 74 by (2.25). 
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Expansions valid in regions complementary to S(0) are (5.11), (5.28), (5.29), (5.30), 
(5.28*), (5.29*), and (5.30*). When O<0<$?r the expansions (5.28) and (5.29) between them 
cover the principal curve ER of figure 2; expansions (5.28*) and (5.30*) cover the principal 
curve AQ, except whenO<0<e and Re£> — 1— 5'. The region of validity of (5.11) however 
includes that part of AQ to the right of Re t> — 1 + 5' when ()<6<e. Hence the expansions 
we have obtained cover the whole £-plane (except of course the neighborhoods of the points 
2=±1) when O<0<j7r, and similarly also when — ^7r<0<O. 

The branches of the functions f, (t 2 — 1)* and ,c/,(£) in (5.28), (5.29), (5.30), and the con- 
jugate results are the same as in (4.3), and the function </(/x) again lias the asymptotic expansions 
(6.1) and (6.2). 

In (5.11) </(ju) is identical with </(/*), and rf=v(t) is given by (5.9). The branch of the latter 
function is continuous in the region of validity of (5.11) and positive when — 1<^<1; the same 
is true of the branches of (1— 1 2 )± and (1— t 2 )%. 

The expansions for the derivative U' {—\\i 2 , pt-^2) corresponding to (4.3) and (5.11) are 
(4.13) and (5.23) respectively, the coefficients &*(%) and v s (t) being given by (4.14) to (4.16). 
The derivative expansions corresponding to (5.28) to (5.30*) have not been recorded, but they 
can be deduced immediately by use of the formal identities 

^7?^£(=F)«=^-±^(^l)*^£(T)'^, (6-3) 

at (t 2 —l)? s=o m *=o M 

and are valid under the same conditions as the parent formulas. 

Uniform asymptotic expansions of solutions of eq (1.1) for large |/*| and unbounded t 
have previously been given by Watson [12], Schwid [10], and Darwin [1]. 

Watson applies the method of steepest descents to integral representations of 



D n (z)mu(-n~>*) 



and obtains a set of asymptotic expansions for the function D n (2^/n) when \n\ is large and 
|arg n\ <J -k— e, which cover the whole f -plane other than the neighborhoods of the points 
f=±l. For each expansion the f -region of validity is fixed, unlike the regions obtained 
in the present paper which are permitted to depend on arg n. For example, the region of 
validity of Watson's expansion analogous to our principal expansion (4.3) is 



UoQ^uu.^t) (|f-i !>«•). 



This is, as we would expect, the region common to S(arg /x) as arg/z varies over the interval 
(— iT+e, \ir— e). 

Schwid applies the methods of Langer for the asymptotic solution of differential equations 
and his results resemble those we have obtained. The main differences are that Schwid confines 
attention to the regions Re 2>0 and |arg /x|<^f7r, gives only the leading terms in each expansion 
and identifies the asymptotic solutions by considering their behavior at t=0 rather than £=oo. 

Schwid 's results apply to the even and odd solutions of the differential equation. An 
effect of carrying out the identification at £=0 is that the error term in the asymptotic repre- 
sentation of each of these solutions is obtained essentially in the form 



e 



- M, *0(|M|- 2 )+^0(|/i|- 2 ) ( 6 ' 4 ) 



in the present notation. If we try to deduce from these results the asymptotic form of the 
function U (— |m 2 , v>t^2) by means of connection formulas, the error term will necessarily also 
be of the form (6.4). Accordingly, a nugatory result will be obtained in the region (U (arg /x)) 
in which U (— i/* 2 , j^t-y/2) is exponentially small for large \,u\. Indeed, the proof given on page 
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358 of Schwid's paper of the fact that the error term in the asymptotic expression for the Her- 
mite polynomials is of the form 6~ M *0(|/i|~ 2 ) instead of (6.4) is false. The mistake is apparently 
due to overlooking the fact that the O's are functions of the independent variable as well as of 
the parameter p. 

Darwin develops series solutions of the eq (2.10), and in [4] Miller applies the same methods 
to eq (2.1). The results are typified by the series 

<7(a,z)= {r( jTf )}i exp S-e-±lnX+f; (_)■ ^\ (6.5) 

(2ir)* L * «=i A J 

for negative a and real positive 2, in which 

X=j(z 2 +±a), d=\ f Xdz, 

L t/2V(-o) 

and d Zs is a polynomial in z of degree 3s if s is odd, or s if s is even. Explicit expressions for 
d Zj d 6 , . . . , d 2 ± are given on page 68 of [4]. 

The series (6.5) is derived by purely formal methods from the differential equation with- 
out any investigation of its asymptotic nature. Applying the substitutions a= — |ju 2 and 
z=nt^2 we find that X=^(t 2 — 1)V2 and 0=/x 2 £(£), in the notation of sections 3 and 4. Next, 
replacing {r(|— a)}i = {r(J-f-|ju 2 )}^ by its asymptotic expansion for large positive n (cf. 
(2.22)) and expanding exp {^Z{—) s d Zs X~* s } in descending powers of X 3 , we find that the series 
(6.5) reduces to the form (4.3). Hence (6.5) must be a uniform asymptotic expansion for 
large \a\ of the same character as (4.3). 

Darwin remarks that the explicit expressions for the coefficients d 3s are simpler in form 
and easier to calculate than the expressions for the coefficients in series of the form (4.3). 
This is indeed true, but we cannot however regard the form (6.5) as being entirely superior 
to (4.3) ; it does not, for example, lead to convenient forms of expansion when we differentiate 
with respect to z. 

Part 3. Expansions in Terms of Airy Functions 

7. Asymptotic Solutions of the Differential Equation 

In this section we seek to apply theorem B of [9] to the differential eq (1.1). The prelimi- 
nary step is to take new variables f=f(£) and W, defined by 

'GO-". H&H^ 

(cf. [7], sec. 2). Equation (1.1) then becomes 



where 



on reduction (cf. (3.3)). 

The first of eqs (7.1) may be integrated to give 



d 2 W 




d 'V**i 5 1 r^m- 


5 (3i 2 +2) 


di) J i6f 2+r>(?j - 


16f 2 4(< 2 -l) 3i 



(7.2) 
(7.3) 



I r*=JV-l)^=|*(* 2 -l)*^ln{*+(* 2 -l)*}, (7.4) 

on specifying that f=0 when t=\. In terms of the variable £ of section 3, 

2 



3 J*=fc (7.5) 
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An essential difference between the functions £(£) and f (t) is that the former lias a branch 
point at t = l, whereas the latter is analytic at this point. Both functions have branch points 
at t = — l. The branch of £(t) which will be the most convenient for our purposes depends on 
0==arg /x. We define it to be positive when t^>l and continuous in the ^-plane cut along the 
principal curve ER of figure 2 when O<0<7r, along the curve conjugate to ER when— ir<0<O, 
or along EQ and ER when 0=0. This determines f(£) uniquely everywhere except in the 
region (U 2 (0)) to the left of EQ+ER when 0=0. 

Following [9], section 1, we make the further substitutions 



0=arg M , t=e~i ie z. 



Equation (7.2) then becomes 



dz 



-— wz- 



-V $ F(e-* i9 z)}W. 



(7.6) 



(7.7) 



The mapping of the £-plane on the 2-plane is conveniently carried out by using x(t) (defined 
in sec. 3) as an intermediate variable; in terms of z 



2 * 

X=- 2*. 
o 



(7.8) 



In our first application of theorem B to eq (7.7) we make a cut in the s-plane along the 
part of the level curve of the function exp (f z%) which joins the point E of affix— (%ir)*e& and 
corresponding to t= — l, to the point at infinity on the negative real 2-axis. Such a join can 
be made when |0| <%w (cf. [7], fig. 2). 




a,(fl) 



Figure 6. z-plane: domains 
D(0) andG(6) (|0|<^r-eY 



Figure 6 illustrates the mapping, the lettering corresponding to figures 2 and 3. The cut 
is the heavy curve EQ. The principal curves through A are mapped on the rays arg z= ±^t 
and 7r; the principal curves through E are mapped on the two sides of the cut EQ and the 
remaining part of the level curve of exp (fs^) passing through E. Cut in this way, the 2-plane 
is mapped on the complement of the domain U 2 (0) (see fig. 4) in the £-plane. 

When |0| < Jt— e, we take the ^-domain D(0) of [9], section 9, to be the cut plane with the 
circle center E and radius J5 removed, 5 being the fixed small positive number introduced in 
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section 3. Then F(e & d z) is regular in D(0) (including 2=0), and the condition (9.2) of [9] 
is satisfied with <?=% because 

e -V»F(e-i te z) ^ (7.9) 

as |s|— >oo 7 uniformly with respect to bounded arg z and bounded 0; this result follows from (7.3), 
(7.5), and (3.5). 

We take the ^-domain G(0) required by [9], section 9, to be the complement of the map of 
the ^-domain V 2 (0), defined in section 5. The boundaries of G(0) are the dotted curves sur- 
rounding the cut EQ in figure 6. All the preliminary conditions of [9], section 9, are then 
satisfied if 5= 5(e) is taken to be small enough. 

Taking the point &i(0) of [9], section 9 to be at infinity on the positive real axis, we find 
Hi(0) = G(0). Hence from theorem B, a solution W\ of (7.7) exists, such that for large positive 
values of u 

w \w % N ^ A s (0,z) . Ai'(uh) • B 8 (B,z) / ^, /M , Q ,A \ 
W 1 ^Ai(u*z)^^t J -+ ^2—727- (seG(0), |0|<^7r-e) 

uniformly with respect to z and 0, the coefficients A s (6, z) and B s (d,z) being given by recurrence 
relations derived from eqs (9.12) and (9.13) of [9], When the variables /x and f are restored, 
this result becomes 

Wt-Aitffit^+^t*® (7.10) 

8=0 M /jj s=0 M 

as |/x|->oo ? when 

t e <7V 2 (0) and -~ 7T+e <6<~ 7T-e, (7.11) 

where A (^) — constant and 

B.(r)^r*fV*{-^(r)+md.(r)}<*r] 

J ° V (8=0,1,2,...), (7.12) 

A, +1 (0=-lB',«)+lj ^(f)B.(f)df J 

the paths of integration lying in the map of CU 2 (0). 

A second application of theorem B can be made to the differential eq (7.7) by taking the 
cut in the 2-plane for the function F(e~^ ld z) to be the part of the level curve of exp (f z%) joining 
E to S, the point at infinity on arg z= — \ir. This join is possible if O<0<7r; the corresponding 
^-domain is then CTJ 3 (0). 

With similar analysis, again taking a x (d) to be the point z=-\-™, we may show that a 
solution of eq (7.2) exists having the right-hand side of (7.10) as its asymptotic expansion for 
large |/x|, uniformly valid with respect to t and when 

t e <7V 3 (0) and e<8<T-e. (7.13) 

The coefficients A s (£) and B s (£) are again given by (7.12), the paths of integration now lying 
in the map of CU 3 (0). 

By considering the limiting behavior as f— >oog~i 1 ^ we see that when e<0< %w— e, the ratio 
of the new solution to W\ is independent of f and has an asymptotic expansion for large |/x| 
identically equal to unity. Hence the expansion (7.10) for W x holds with the conditions (7.13). 
Similar analysis, or an appeal to symmetry, shows that this expansion also holds with the further 
conditions 

t e C r V 1 (0) and -7r + e <0<-e. (7.14) 
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We can combine the three ^-regions of validity into a single region T(0), given by 

CVi, CV^CYo, C\ 2 , CY 2 \JCY Z , <7V 3 , 

according as d lies in the respective intervals 



(— 7T+e, — £ 7r+e), (— 2 **+€, 



-€)j ( — €, e), (e, ^tt— e), (^7T— €,7T — e). 



(7.15) 



T(0) is illustrated in figure 7 for O<0<7r— e; T(— 0) is the conjugate of T(0). The dotted curves 
are of the same character as the dotted curves of figure 2; their z-maps parallel the z-maps 
of the principal curves through E at a distance 8 away. T(0) is the unshaded domain bounded 
by the dotted curves. 

Collecting together the results of this section, we have shown that there exists a solution 
Wi of the differential eq (7.2) whose asymptotic expansion for large |/x| is given by (7.10) and 
is uniformly valid with respect to t e T(0) and in the interval (— 7r+e,7r— e). The functions 
£(t),F([), A 8 (?) and 2?,(f) are given by (7.4), (7.3), and (7.12), and are regular when t e T(0). 




(a) 0< 6 < e 



(b) €<9< fir-e 



(c)-^7r-e< 6< it -e 



Figure 7. t-plane: domains T(0). 
8. Identification of Solutions 
From (4.1) and the second of eqs (7.1), we see that 



(8.1) 



satisfies the differential eq (7.2). When \x is fixed and t-+ <» e ~ , ' arg " this solution becomes ex- 
ponentially small. If \n\ is sufficiently large and |arg /x| <w— e, the solution M^of section 7 is 
also exponentially small in the same circumstances, for from (7.4) we have f~(f) T ^, and 

4 

so/x^f-^+°°. Hence the two solutions are linearly dependent. Accordingly, 

u(-±^)~m (£f[* o* g^M^g&jp] ( , 2) 

as |/x|— >oo, uniformly with respect to t e T(arg /x) and arg/x in the interval (— 7r+€, 7r— -e). 

This is the fundamental expansion in terms of Airy functions. In order to determine the 
function h(n) y we first fix the arbitrary constants of integration in the recurrence relations 
(7.12) by specifying the conditions 



A(D = 1, lim4,(f)=0 («=1,2, 
iri->« 
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)• 



(8.3) 



We now examine the limiting form of the ratio of the two sides of the relation (8.2) as t-> <» e~ i{iTg,i - 
With the aid of (2.2), (3.5), (7.5), (8.3), and the asymptotic forms 

Ai &**f)~i *"*(f**f)" l «" rt , Ai'(M*f) 1 x-J^f)*^* 

(cf. [8], appendix), we find that l/A(/x) has the asymptotic expansion 

as |/x| — >oo ) uniformly valid with respect to arg /x, where 



(8.4) 



A.= ]im{r*£,(f)}. (8-5) 

lfl->« 

An alternative form of expansion, not involving the quantities h s , may be found as follows. 
Suppose that e<0< Jir— e. The principal curve AQ is then included in the domain T(0). 

4 

On this curve /* T f is negative (cf. fig. 6). Let us set 

r=f^/T*<4 (8.6) 

where gC is the mth negative zero of the function Ai'. Then (8.2) becomes 

U (-i m 2 , mW2)~A (m) (jTTi)* Ai (O g ^£H (8.7) 

in which £ m denotes the value of t corresponding to f=f TO . Now keep ju fixed and let m->oo, 
so that t-><x>e . Then from [8], appendix, we have 

<£=-{fir(4m-3)y {l + 0(m- 2 )}, Ai (0= (~) w " V~ § (-(4)"* {l+0(m" 2 ) }, (8.8) 

and from (7.5) and (3.5) 

&-ln2* m =| £+|+0(| * m |~ 2 ). (8.9) 

Substituting in (2.7) by means of these equations, we find that the limiting form of the left- 
hand side of (8.7) is minus the expression (5.17), with t m defined of course as in the present sec- 
tion. On the other hand, the limiting form of the right of (8.7) is 



-*- ' 



(cf. (8.3) and (8.8)). Comparing the two sides we find that \ 7r~i /x"* h(n) has the same asymp- 
totic expansion for large |/x| as the right-hand side of the relation (5.18). Hence we obtain 

A( m )^27tM.?(m)^2^M#G0, (8.10) 

where gin), g(n) are defined in sections 4 and 5; g{y) has the asymptotic expansions (6.1) and 
(6.2) for large |ju|. Although (8.10) has been proved on the assumption that e<arg id<%Tr—e, 
(8.4) shows that l/A(/x), and hence also A(/z), possesses a single asymptotic expansion over the 
range |arg ju| <w— e. Therefore (8.10) must also hold over this wider range of arg ju. 

Substituting (8.10) in (8.2), we obtain the final form of the fundamental expansion, given by 

V (-1 ,»,< ^ * „ (£)* [Ai („.„ £ ^2+M^C § M>] (8 .u> 

as |/i|— »oo, uniformly with respect to i e T(arg p), arg ju in the interval (— ?r+«, 7r— *). 
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The coefficients A s ({) and B s (£) appearing in (8.11) are determined by the recurrence 
relations (7.12) and the conditions (8.3). Explicit expressions can be found by the method 
used in [8], section 6 for the determination of the coefficients in the uniform expansions of Bessel 
functions of large order. We take f to be fixed and positive, and let /x-^ + co . The Airy func- 
tions appearing in (8.11) can then be expanded in their asymptotic series for large positive 
arguments and the result compared with (4.3). Thus we derive the expressions 

Is 2s-4- 1 

^=23^^(0, ^(D-S^^te (8-i2) 

m=0 ra=0 

(cf. [8], (6.6)), where ja4(£) is defined in section 4, Oo=l and 

„ _ (2m + l) (2m+3) ... (6m-l) . 6m + l „ , filQ v 

am ~ m!(144)~ m 6m-l flm ' { * A6) 

Analytic continuation shows that eqs (8.12) hold for all values of f. 

We note from (8.10), (8.4), and (6.1) that 

A»-lm{rte.(r)}=-Jto+i (s=0,l, . . .), (8.14) 

iri->- 

where g s is given by (4.12) and (5.21). 

Expansion for the derivative. The expansion (8. 11) may be differentiated t erm by term with 
respect to t; this is a consequence of theorem B. With the aid of the first of eqs (7.1), we find 
that 

V 2 / 0(f) [_ m t 5 =o M 5=0 M J 

as |ju|^°o, uniformly with respect to t e T(arg ji) and arg p. in the interval ( — 7r + e,7r— e). 
Here we have used the notation 

*<» -0^)41)' 

The coefficients C,(f) and Z? s (f) are expressed in terms of A„({) and -B s (f) by the equations 

c,(r)=x(fu 8 (f)+^(f)+f5,(f), z? I (f)=it(r)+x(f)B.. 1 (f)+B;_i(f) > (8.17) 

(cf. [8], (6.10)), where 

X(f)= ^y _ 4f-2^ziyi- -*r~" (8 - 8) 

Explicitly, 

r*ft(r)=-S ^r* w ^2.-«+i(*), z?.(f)=s a w r* w ^2.-m(«), ( 8 - 19 ) 

?n=0 ra=0 

(cf. [8], (6.12)), where #*(£) is given by (4.14) to (4.16), and a m , b m by (8.13). 

9. Use of Connection Formulas 

In this section we seek expansions which hold in /-regions complementary to T(arg /*) when 
|arg n\<i ir (cf. sec 5). 

In (8.11) replace /x by i*e ir . Then using (5.24) we obtain 

V (~l ^,-m(V2)~2»'«"*"-' Vi/M (^i)' 

L »-o A* ^"/i* s=0 m J 
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when 

-T+e<d<-€, teT(d+T). (9.2) 

In particular, we see from (7.15) and (5.3) that (9.1) holds when 

-~7T<0<-e, t eCY 3 (d+w) = C\ 1 (d). (9.3) 

Similarly, 

X [Ai (.1 VI) S ^#+^1^ £ ^#1. (9.4) 

L 5 = M g 3 /X 7 * = M J 

when 

€<^<^tt, teCY z (6). (9.5) 

The expansions (9.1) and (9.4) are of course the expansions for U(— J/z 2 , m^V2) in terms 
of Airy functions valid at the turning point t= — 1. Combined with (8.11) they cover the 
whole of the Z-plane when e< |0| <\ w. 

When \6\ <e we cannot expect to be able to represent the function U(—i/x 2 , fxt-y/2) in 
Y 2 (6) = CT(d) by means of a single Airy function, because as \t\ — >oo in this region U (— -J/x 2 , /^V2) 
becomes exponentially small or exponentially large according as /x 2 is or is not an odd positive 
integer (cf. (2.7)). 

To deal with this remaining case we use the connection formula 

17(4* - M * V 2)-siniv-tf(-iM», M^)=@ i g^C/QM 2 ,-^V2) 

+(I) rcFpy^G"'^ 2 ) (9 - 6) 

obtained in an obvious manner from (5.26). Replacing \x by ^e~^ lir in (8.11) and using the first 
of eqs (5.25), we find 

We substitute this result and its conjugate form on the right of (9.6), and on the left we sub- 
stitute (8.11). Using the eq [3] 

Bi ( M *f)=<rHAi (e -+Vf) +*** Ai («*Vf) . 
and its differentiated form, we obtain 

Z ^ S = M jLt 3 5 = M J J 

This is valid when \9\<^w—e in the £ -region common to the domains T(0), T(0— J x), and 
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T(0+^7r). Using (7.15) and (5.2), we see that this means (9.7) is valid when 

|0|<|x-€, teC\ 2 (d). (9.8) 

This is the required result. We note that when \x 2 is an odd integer the terms in Bi, Bi' 
in (9.7) vanish, as we require from a consideration of the behavior of U (— |m 2 > m^V2) as |£|^°° 
in the sector V 2 (0). We may also verify without difficulty that (9.7) agrees with (9.1) and 
(9.4) in their common ^-regions of validity when — ^7r+e< (9< — e and e<6< } 2 tt— e, respectively. 

10. Summary of Expansions in Terms of Airy Functions 

The fundamental expansion of this kind is (8.11). It is valid for large \/jl\ uniformly 
with respect to = arg \x and t, when — 7r+e<0<7r— e and t lies in the domain T(0), defined 
by (7.15) and illustrated in figure 7. The function #(//) is calculable from either of the asymp- 
totic expansions (6.1) and (6.2). The function f (t) is given by (7.4), that branch being chosen 
which is continuous in T(0) and positive when £>1; similarly the branch of {{/(t 2 — 1)}* is con- 
tinuous in T(0) and positive when t^>l. The coefficients A s (t) and B s (£) are given by eqs 
(8.12); they satisfy the recurrence relations (7.12) and the eqs (8.3). They are regular func- 
tions of f in the f-map of T(0) and real when £* lies on the real axis to the right of the point 
of affix -(fir)*. 

Other asymptotic expansions for U(— h^ 2 , a^V2) are (9.1), (9.4), and (9.7); they are 
valid with the conditions (9.3), (9.5), and (9.8), respectively. The combination of these 
expansions and (8.11) covers the whole £-plane when |0|<i7r. 

The expansion for the derivative U'(— -|/x 2 , /x£y2) corresponding to (8.11) is (8.15), and 
holds under the same conditions. The coefficients C s (f) and D s ($) are given by (8.17) and 
(8.19). The derivative expansions corresponding to (9.1), (9.4), and (9.7) are not recorded, 
but they can be obtained immediately by analogy with (8.15). They are valid under the 
same conditions as the parent formulas. 

Approximations for parabolic cylinder functions of large orders which hold in regions 
containing one of the turning points have been given by Watson ([12], sec. 17), Schwid ([10], 
sec. 6), Erdelyi, Kennedy, and McGregor [2], and KazarmofT [15]. 

Watson derives asymptotic expansions in terms of elementary functions for the functions 
D n (z±2^n) = [/(— n— J, z±2^n) as |tx| — >oo ? with the conditions |arg n|<i7r and z=o(\n\~$). 
In the present notation these conditions correspond to |arg ju|< iTrand^Tl^od/xl - ^). Watson's 
expansions are analogous to MeisseFs expansions for Bessel functions of large order and may 
be obtained from the present uniform expansions by suitable re-expansion (cf. [8], sec. 6). 

Schwid gives approximations in terms of Bessel functions of order one-third (Airy func- 
tions) which are valid in the neighborhood t—l = 0(\ii\~ 2 ) as |/x|— >oo, using the present notation. 
This neighborhood is complementary to the region of validity of Sehwid's approximations 
in terms of elementary functions, described in section 6. No investigation is made of the 
validity of the Bessel function approximations outside the neighborhood. 

The approach of Erdelyi, Kennedy, and McGregor is similar to the one we have used. 
These writers obtain a set of asymptotic approximations in terms of Airy functions, to solu- 
tions of eq (1.1) when \x and t are both complex (using the present notation). The results 
are contained in those given in the present paper; the principal ways in which they have been 
extended and improved upon are as follows. 

First, the approximations apply to fixed values of arg /x, no investigation is made of their 
uniform validity with respect to arg /x. 

Secondly, the approximations are established with the condition |arg^|<|7r; outside this 
interval connection formulas are to be used. This leads to complexity in the number of results; 
eleven are given compared with the four we have found to be sufficient. 

Finally, only the leading terms of the asymptotic expansions are found; no terms involving 
the first derivative of the Airy functions appear in the approximations. An effect of this is 
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to exclude regions in the ^-plane containing the zeros of the Airy functions. No special treat- 
ment of such regions is required in the present theory. 

KazarinofTs paper (which has appeared since the present paper was written) contains 
a brief investigation of the extension of the results of Erdelyi, Kennedy, and McGregor to 
include uniformity with respect to arg /x and higher terms in the asymptotic expansions. The 
results are rather more complicated and less complete than those given here. 

Part 4. Results for Real Variables 

In this part we suppose jjl, t, a } and x to denote real variables except where otherwise indi- 
cated; n is alwa}^s large and positive. The variable x is not to be confused with the complex 
function x(t) introduced in parts 2 and 3. 

11. The Equation w" = (^x 2 -\-a) w 

The principal solutions of this equation are U(a,x), U(a,—x), V(a,x), and U{a ) x) J defined 
in section 2. 

(i) a positive. When a — > + oo the differential equation has no turning points on the 
real x-axis. The expansions of greatest interest are therefore those in terms of elementary 
functions, and we shall not consider the forms of the Airy expansions. 

In (4.3) replace n,t by ip, —it, respectively. Then using (4.8), we obtain 

(1 \ 1- z?M 2 £(-*0 oo y (—if) 1 

as fx — > +oo. In the complex plane the region of validity of this expansion is e? i1c S(%T), and 
from figure 2(c) it is clear that this region includes the whole of the real £-axis. 

We have now to interpret the branches of the various functions. 

From (6.1) and (6.2) we derive 

g(^)=e-^^ 2) g(fji), (11.2) 

where the function g(fx) is real and has the expansions 

gW^'-h^-iS-hfl+i^i-y-^X, (11.4) 

I 8=1 (lM ) J 

as ix -> +oo. Comparison of (11.4) with (2.27) shows incidentally that 

l„_l u 2_l _l u 2 l u 2_l 
-/ X r.li.2-4 l u 2 _l n 2_l , T 2 2 4M *e * M M 2 2 m ~\ 

000^2*" fe* M M 2 * 2 H r a-Li *T ' (1L5 ) 

1 V2 I 2M ) 

where the symbol c^z is defined in section 5. 

Next, when t lies on the negative imaginary axis of figure 2(c), £(t) lies on the straight line 
O'C shown in the third diagram of figure 1. Hence 

f(— iO=|iir— 1(0, (H.6) 

where J(t) is positive when £>0. From (3.4) we find 

~ m= Jo ^ +l ^ dt= \ t(* 2 +l)*+| ln {*+(* 2 +l)l}, (H.7) 

the branches here having their principal values. 
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Lastly, we have 

(-t 2 -l)i=e-i l7r (t 2 +l)*, (-t 2 -l)% s =e-% si "(t 2 +l)$ s , (11.8) 

where the functions on the right have their principal values. 

Substituting in (11.1) by means of (11.2), (11.6), (11.8), and the relation 

u s (t)^i s u s (-it), (11.9) 

we obtain the required result in real form, given by 

(j(LvV2V^W- ^rS Us(t \ 4 (- 00 <K °), (11-10) 

uniformly with respect to t. Explicit expressions for the coefficients u s (t), obtained from (4.11), 
are 

u Q (t) = l, u 1 (t) = (-t 3 -6t)l24, iZ 2 (£)== (9* 4 + 249/ 2 - 145) /l 152,1 

y (li.ii) 

12 3 (0 = (4042£ 9 +18189Z 7 +28287£ 5 -l 51995£ 3 + 2 59290^/4 14720. J 
The corresponding result for the derivative is 

valid with the same conditions. Here 

v s (t)=i*v 8 (-it). (11.13) 

Explicitly, 

^(«) = 1, 5i(0=(— ^+6«)/24, ^ 2 (0 = (— 15^ 4 — 327^ 2 +143)/1152, 
£ 3 (/) = (4042£ 9 +18189£ 7 +36387£ 5 +2 38425£ 3 -2 5929(M)/4 14720, 

(cf. (4.13) to (4.16)). 

Since U(a,x) and U(a, — x) comprise a satisfactory pair of solutions of the differential 
equation when a is positive, it is unnecessary to give the corresponding expansions for the 
solutions V(a,x) and U(a,x). 

(ii) a negative; expansions in terms of elementary functions. For t>l + 8, the expansion of 
£7(— Jm 2 ? /A/2) is given in real form by (4.3 ) 3 . In this expansion q(jjl) is calculable from (6.1), 
(6.2), or the relation 

g(fi)a- 2-^ 2 -V^V iM2 ~H^2^ 2 -*^ (11.15) 

and £, j^(£), by use of (3.4), (4.8), and (4.11), all fractional powers now taking their principal 
values. 

When — l+5<t<l — 5, the expansion of U(— |/x 2 , ^V2) is given by (5.11). Here g(v) is 
identical with gin) and rj is given by (5.9). 

For the remaining part of the real axis we have the expansion (5.28). In the present 
circumstances the function £(t) is positive and bounded away from zero. Hence the contribu- 
tion of the series i cos hnx 2 • i~ M * y^, s£(£)u~ 2s is always asymptotically negligible compared 
with that of 2 cos V*^I](-) s 4(Qm"" 2s . Thus (5.28) reduces to 

t/(-Lv-W2)~ -M-\ 2 cosi *fi. + ± i-)>^+ S ml>.e-*±,^m (11.16-) 

\ z / it — 1)^ L L s=0 M * *=0 M J 

3 For convenience the symbol 5' of parts 2 and 3, denoting an arbitrary fixed positive number, is replaced here by 8. 
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(11.14) 



valid when t>l + 8. The functions g(n), £ andj^(£) are the same as in (4.3), discussed above. 
The corresponding expansions for the second solution U(a, x) are of interest, particularly 
when x>0. It is not possible to obtain them trivially by substituting in the connection formula 
(2.12) the expansions we have just obtained for U(a, x) and U(a,—x), because although can- 
cellation of terms in the asymptotic series takes place, we cannot assert immediately that the 
remainder terms implied in these series also cancel. To put the results on a sound basis we 
substitute the expansions (5.5) and (5.6) in the formula 

u (-§ ,\ W2)= (2.) -* r Q+i /) [r-a+w u Q ^-^2) W^ +i * 2) u Q ^ ^2)} 

(11.17) 
obtained from (2.5) and (2.12). 

When £>l + 5 the branches of £(t) and jz/ s (£) m both expansions are the principal ones. 
With the aid of the relation 

(2t)-* r Q+^m 2 ) e^i+W gfa^cxgQi), (11.I8) 

proved by use of (2.22), (5.7), (5.20), and (6.2), we deduce that 

Tj(-lix 2 ,W2V 2 ff( m)-^—S (~) s M ' ffl 3 4 (*>l+«). (11.19) 

V 2^' MV / * vw (J 2 — l)*£o y ($*—l)* a M* - ^ ) } 

When —l-\-d<t<l — 8, we interpret the branches of the many- valued functions by use 
of (5.8) and (5.10). With the aid of (11.18) we obtain on reduction 

U[ — -/x , /xrV2 J^ r cos(m r?4-x7r) > , — j9x „, , ~sm (/oi-t-z7r) > f -^ > 

V 2^ ,M / (!—$»)* L ; ^(l-£ 2 ) 3 V 4s ^o(l-^ 2 ) 3s +V s+2 J 

(11.20) 
valid when — l + <5<£<l-<5. 

The expansions for the derivatives corresponding to the expansions obtained in this sub- 
section can be written down immediately by use of (6.3), or by analogy with (5.23). 

(iii) a negative; expansions in terms of Airy functions. For t> — 1 + 5 the expansions of 
U(— -iV,^V2) and U{—^ix 2 —^2) are given by (8.11) and (9.7), respectively. In these 
expansions f (t) is given by (7.4) and is real when t^> — 1 ; a more convenient form of this relation 
when — 1<^<0 is 



c— r)*=J] 1 (i— o*rf*=| cos-**-- *(i-*)*=ij(o (n.21) 



(cf. (5.9)). The coefficients A s (t) and B s (f) are real and given by (8.12), (8.13), and (4.8) 4 . 

An expansion for the solution U(a,x) may be derived by use of (11.17), as in the case of 
the exponential-type expansions. The result is 

^(-^^^V2)~2.V,w(/ T )TBi(M^)S^+^¥^2^#} (11-22) 

valid when t > — 1 + 5. 

The corresponding expansion for U'(—iiJ?,nt^2) may be written down by analogy with 
(8.15). 

(iv) Hermite polynomials of large degree. Expansions for the Hermite polynomial H n (x) 



* More convenient expressions when — 1<£<1 are given by (13.4) and (13.5). 
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when n is large can be obtained from the results in subsections (ii) and (iii) above by use of 
the equation 

H h ^t)=2W-^ n% u(-\ mVV2} (11.23) 

Thus (4.3) and (5.11) yield expansions which hold uniformly in the ranges #>(l + 5) (2// + 1) 2 
and —(1 — 8) (2n-\-l)* <x < (1 — 5) (271+1)*, respectively. An expansion uniformly valid through- 
out the range x> — (1—5) (2n+lp is furnished by the Airy expansion (8.11). 

12. The Equation w"=(a—\x 2 )w 

(i) a positive; expansions in terms of elementary junctions. From (2.16), (2.17) and (2.18), 
we have 



where 



k~*wQ ^ 2 ,^2)+i0wQ /x 2 ,- M fV2)-2M^V ( ^+i-)[/ Q ifo-frtitTpy (12.1) 

k=(l + e^-e^ 2 , ^' 02 = ±{ r Q+^'^ ( 12 - 2 ) 

For large ji, 

k= ^ e -l^ + 0( e -^ 2 ), ^^=6-^ 2 Qm 2 V V {1 + 0(m- 2 )}- (12.3) 

In (4.3) replace /jl by fie ~^ 1 . Then we obtain 

U (I i^e-^my,^-) ^L, S «■ ^ jr. (12.4) 

as M"^+ oo . The £ -region of validity of this expansion, S(— \t), includes the segments t> 1 + 6 
and -l + 6<^<l-5 of the real axis (cf. fig. 2(b)). 
When /i^+co, we have from (12.2) and (2.22) 



e* l + 2 ~e 



-w A „ 2 \ jiM 2 f^ y* 1 *| ^ 7» 
V2 M / L£<S(W)'J |« (W 



(12.5) 



the choice of ambiguous sign on taking the square root being resolved by use of (12.3). From 
(6.1) 

giiur**) ~2^ 2 -iel i "V''-i"'V i4 ' 2 ~ i (l+ / ii)" 1 , (12.6) 

where the symbol L denotes the asymptotic series 

I-Z}(-) ,S £& (12.7) 

Further from (5.21), with n replaced by /le 1 ", we have 

,ttf(fiM 2 ) 8 l-i£ 1-iL' l ; 

Introducing a function l{n), defined by 

l(n) ^V^V^-^Gue -1 * 1 ), (12.9) 

and using (12.5), (12.6) and (12.8), we find that 

., , 2i/l±iL\i\l+iL-i 1 2* 1 f . 

l( ^ ~ -AT=Tl) \i=iL 1+3 ~ ? a+Zr)V (12 - 10) 
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Thus 1(h) is a real function, having an asymptotic expansion of the form 



?(j")~^i:4 m m-*+»- (i2.il) 

{JL 2 s = M 

From (12.10), (12.7), and (4.12), we find 

We are now in a position to combine (12.1) with (12.4). When t>l + 5 each of the func- 
tions £(£), (t 2 —l)* and (£ 2 — l)^ s is positive. Substituting by means of (12.9) and the first of 
(12.3), and equating real and imaginary parts, we find 

urn 2 s k>\ 2~V* 7r/x2 /(M) T ( 2i.li \ ^ (-Y u 2s(t) . M ,, n A (-)% 2s+1 (0 "1 

(12.13) 

w/i 2 , /o\ sV^r. , , ^ (-) s u 2s (t) , , 2 . , _ , ^ (-)%,+iW 1 

w(iM 2 ,-^v2) — ^=i)F L sm (M m7r) § (^i)v +cos (M m?r) S ( ^i)3^V +2 } 

(12.14) 

uniformly with respect to t when t > 1 + 5. The function Z(/x) is calculable from (12.11). 

When — 1 + 5< / <1 — 8 the appropriate branches of £ and (£ 2 — 1)* are given by %=—ir} and 
(£ 2 — l)*=e Ixi (l— £ 2 )±, where 77 is the real function (5.9) and (1— Z 2 ) 1 is positive. Hence we derive 

W(L?,M2)~ fy, **' ^^"f i (— 1+*<*<1— *). (12.15) 

\2 ' / 2^ 7r/l2 (l-^ 2 )^1^o (1— p)* f /* ff 

The expansions for the derivatives corresponding to (12.13), (12.14), and (12.15) may be 
obtained by term-by-term differentiation with respect to t (cf. (5.23) and (4.13)). 

It may be noted that (12.13), (12.14), and (12.15) are respectively equivalent to the formal 
series (317), (318) and (324) given on pages 84 and 85 of [4] (cf. sec. 6). 

(ii) a positive; expansions in terms of Airy functions . In (8.11) replace ix by ^e~ Jirt . Then 
we obtain 

The i-region of validity of this expansion, T(— Jx), includes the segment t>— 1 + 5 (cf. fig. 
7(b)). Now from [3] 

2«*-*Ai (e-*Vf)=Bi (- M *f)+iAi (- M *f). (12.17) 

Hence substituting (12.16) in (12.1) and equating real parts, we find, with the aid of (12.9) 
and the first of (12.3), 

pp^,^--^^-^ Lft(- M Tfl§H -^+ — -j g(-) -^sr-J (12.18) 

valid when t> — 1+5. The corresponding expansion for the derivative is 
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xr-5it=^£(-)-^+Bi'(-,«r)a(-)'^#l (i2.i9) 

(cf. (8.15)). 

The expansion for the second solution may be obtained by equating imaginary parts. 
This gives 

(12.20) 

However, this result is only proved in this way for t>l, i.e., non-negative values of f. This 
is because the error term on truncating the series inside the square brackets at the rath terms is 

When f is negative the terms in Bi and Bi' cannot be absorbed in the terms in Ai and Ai' be- 
cause they are exponentially larger. 

That (12.20) is' valid in the wider range t>— 1 + 5 can be proved as follows. In (8.11) 
replace m by pe?* 1 . Then we obtain 



U Q if, e *•* ,W2)~2,rM*VW*' i ) (pzr)* 

xTAi(-M»f) S (-)« ^#+^=^ £ (-)« W\ (12-21) 



valid when t e T(f 7r)and so including £> — 1 + 6 (cf. fig. 7). Equation (12.1), with t replaced 
by —t, shows that the left-hand side equals 

2-V+"V*(l*i+W [V^(J n 2 -nt-yj2)+ikiwQ ju 2 , W^)l • (12.22) 

Also from (6.1) and (12.7) we derive 

^i^-i^G^i**) 

=2-M* M V'(** a "HMZ(/i) (12.23) 

(cf. (12.6) and (12.9)). The expansion (12.20) may now be obtained by substituting (12.22) 
and (12.23) in (12.21), equating real parts and using the first of (12.3). 

(iii) a negative. When a is negative the differential equation has no turning points on the 
real axis and we may confine our attention to the exponential-type expansions. 

From (2.16), (2.17) and (2.18), we have 

k~iw(-l m 2 , /tfVaVa* w (-\ m 2 ,-^V2)=2*e"i^V(^+w u(-\ ^v^Wa)* (12.24) 

where 

k={l+e-^-e-^ 2 =l + 0{e-^ 2 ), (12.25) 
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and 



,**-± / r (i-| **»)/ I r Q-i iM 2 )| y=e& 2 Qm 2 )"^ 2 { 1 + 0(m" 2 ) } (12.26) 



(cf. (12.2)and(12.3)). 

In (4.3) replace /z, t by m^ 1 ^, te~" 5lrt , respectively. Then we obtain 

<-4^"^^)^0-»-) l-^-^gC-*)-^-^^ (12.27) 

The region of validity of this expansion, ^^S^tt), includes the whole of the real £-axis; this is 
clear from figures 2(b), 3 (hi) and 1. 

Substituting (12.25), (12.26), and (12.27) in (12.24) in the manner of (i) and (ii) and 
determining the branches of the many- valued functions, we obtain finally 



W (- J M 2 ,M<V2) ~-M- [cos Gu 2 I+ fcr) T, 



(-) s % s «) 



« 2 +l)iL v ^' *" y ^ott 2 +l)V* 

^^§ (^' } (12 - 28) 

uniformly when — oo<^<co. Here l=l(t) and u s (t) are given by (11.7), (11.9), and (11.11), 
and Z(/x) is calculable from (12.11). 

The corresponding expansion for the derivative is 



(-) s v 2s (t) 



W (-iM 2 ,^V2) ^( M )(*2 + 1) ^ sin ( M 2| + I,,) g _^ 

+ COS ( M 2 i+i7r) 1] (-)^2,+i(0 1 (1229) 



# s (£) being given by (11.13) and (11.14). 

13. Numerical Examples 

Here we illustrate the use of the Airy-type expansions by two examples. It may be 
remarked that applications would be greatly facilitated by the preparation of suitable numerical 
tables of the functions jr( M ), IQi), f(«), {f/(^— 1)}* and 4,(£), 5,(f), <7 5 (f), Z?,(f) (*=0,1, . . .). 
In the present examples these functions are calculated directly from formulas given in this 
paper. 

It will appear from the examples that the early coefficients in the expansions decrease 
in magnitude as ,9 increases. This property, shared with uniform asymptotic expansions for 
Bessel functions of large order ([8], sec. 6), makes the expansions well suited to numerical work. 

Example 1, To evaluate W(a,±x), W'(a,±x) for a =8, x=7. 

We use the expansions (12.18), (12.19), aud (12.20). Here 

M =4, *=7/(4V2) = 1.23743 6867. 
Hence 

2*6^^=4.05527 596X10 5 , (13.1) 

and from (12.11), (12.12) 

Z(/*)~(2 l /2) (1-0.00000 3391-0.00000 0006+ . . .) = 0.59460 1538, 
so that 

ttVZ(m) = 1.67296 798. (13.2) 

Next, from (7.4) we have 

ft=|^ 2 -l)*-|ln{i+a 2 -l)^}=0.16932 9581. 
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Hence 

f=0.30607 0793, f*=0.55323 6652, { f/(^ 2 — 1 ) } i =0.87122 5817. (13.3) 

From (4.8) and (4.11) we compute 

j*$=l, j^= -0.59504 451, j/ 2 = 2.92479 70, j*J= -26.48264 7, 

and from (8.13) 

_ 1 _5 _385 _ 85085 

oo-l, a x —jg fl-2-4gQg^ a3 ~6 63552 7 

. / _ 7 a i^ a 95095 

Oo-l, *i— jg' ^—4608' ° 3 ~~6 63552* 

Substitution in (8.12) yields 

4>(r) = l, ^i(f) = -0.00649 0, S (f) = -0.03637 97, ^(f) =0.01202. 

Next, we have /**f= 1.94342 840, and by interpolation in table VII of [3] 

Bi(— M*f) = -0.39524 090, Bi'(-M*f)= 0.32386 134, 

Ai(— /i*f)=0.26162 602, Ai / (-/x*f) = 0.59097 966. 

The content of the square brackets on the right of (12.18) is accordingly 

-0.39524 090(1+0.00002 535+ . . .) + 4"*X0.32386 134(— 0.03637 97-0.00004 70+ . . .) 

-0.39554 353. 

Substituting this value and (13.1), (13.2), (13.3) in (12.18), we obtain 

w 1.67296 798X0.87122 5817X0.39554 353 _ ftU21fi4WYl(l -5 
14(8,/) 4.05527 596 X10 5 0.14216485X10 . 

Similarly, 

TF(8,— 7) = 1.67296 798X0.87122 5817X0.26109 871 X4.05527 596X10 6 = 1. 54327 57X10 5 . 
For the derivatives, we have first 

2-M/A(m) =1.87784 306. 

From (4.14) and (4.16), we compute 

#o=l, #1 = 1.00283 725, # 2 = -3.52326 59, # 3 =29.91942 3, 

and thence from (8.19) 

Co(f)= -0.07833 70, Ci(f) =0.00971, D (f) = l, A (f) =0.00760 4. 

Hence 

tit/zotn 1.87784 306X0.31897 315 A .^.n « 19vtA _5 

TF(8 ' 7)= - 4.05527 596X10 5 X0.87122 5817 == -°- 16953613X10 ' 

d"(8 7) _ l-87784 306X0.59419 14.3X4.05527 596X10 5 __ 5 19368 20xl0 5 

0.87122 5817 

The results of this example agree adequately with the entries given in [4], table I. 
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Example 2. To evaluate the Hermite polynomial H n (x) for n=10, x=i by means of the 
expansion (8.11). 

Here (cf. (11.23)) 

M=V21 =4.58257 569, t=4/^21 =0.87287 1561 . 

From (6.1) and (4.12) 

1/^(m)-2 5 -V- 25 (21)- 5 (1+0.00198 4127-0.00000 1052+ . . . ) = 2. 11580 149X10~ 3 , 

and hence 

2*-Mff(/0 =2.78291 471 X10 3 . 



Next, from (11.21) 



A_ 3 __ tJ 3 



(-f) 2 =^ cos-^«-^(l-^) ^=0.06286 6477. 



Hence 

f=-0.15810 5174, (-f)*=0.39762 4413, {f/(£ 2 -l) }*=0.90271 123. 



When f is negative, real forms of equations (8.12) and (4.8) are 



i,(fl = (-)'Si(-f)" !M 4,-., (-f) 5 -B s (f) = (-)'- 1 Sa m (-f) *"jafc_ + „ (13.4) 
where 

^f s =u s (t)/(l-t 2 y s ' (13.5) 

In the present example we find 

J5;=l, ^"=-1.63978 319, ^ = 21.1903125, J^= -516.52363, 
and thence 

-4o(r)=l, A(f) = -0.01024 3, 5 (r) = -0.04317 57, J5,(f) =0.0166. 

Next, /x*f= — 1.20344 324. Interpolation in table I of [3], yields 

Ai0**f)=0.52582 209, Ai'(/i*f) = 0.10920 810. 
Substituting these results in (8.11), we obtain 



I7(— T, 4V2)=2.78291 471X10 3 X0.90271 123X0.52572 856 = 1.32071 865X10 3 , 
and thence using (11.23) 

H l0 (4) = 2*<?u(^-~' 4V2^= 1.25984 22X 10 8 . 

The exact value of this polynomial, computed from explicit expressions given on page 75 
of [4], is 1259 84224. 

Part 5. Zeros and Associated Values 

14. Zeros of U(a, z), 17' (a, z), U(a, z), and I?' (a, z) 

(i) a negative;— 2(1 — 8)^J(—a)<z<2-y/(— a). Let the sth real zeros, counted in de- 
scending order away from the point z=2^/(—a), of the functions U(a,z) } U'(a,z), U(a,z), and 
U'(a,z) be denoted by u a , s , u' atS , u atS} and u atSj respectively. 

The expansion (8.11) is of exactly the same form as the expansion (4.24) of [8]. The 
analysis of section 7 of that paper is therefore immediately applicable, and we find 

« a , s ~2i M L(«)+^+^+...l (14.D 

L MM J 
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as /x=V(~~2a)->+ oo, s fixed. Here 

a=fjL-*a St (14.2) 

a s denoting the sth negative zero of the function Ai and not to be confused with the a m of (8.13). 
The coefficients p r (f) are given by 

».(«=*(«, ft(«=«i(f)*'(f), ^(r)=«2(f)t'(r)+|{« I (f)} 2 r(f), ..., (i4.3) 

where f(f) is the function inverse to f(0, defined by (11.21), and the a r (£) are given by rela- 
tions of the form of (7.13) of [8]. Explicitly, 

^ = 24(*'-l)' + 48 {(*2_ 1){ « } l- ( 14 - 4 ^ 

The corresponding expansion for the associated value U'(a,u a<s ) can be found by sub- 
stituting (14.1) in (8.15). More convenient expressions for the coefficients are obtained how- 
ever, if we substitute instead in the equation 

U'(a,u a ,,) = ±(2*)l{T Q-o)}* (-^f )~* ( 14 - 5 ) 

obtained from eq (2.15) of [6] and the second of (2.13) in the present paper. This yields the 
required result in the form 

U'(a,u a , s ) ~r* M * | r Q+i m 2 )} * Ai'(a,)^(«) [l+ZJ ^} (14.6) 

where 

and the coefficients P r (f ) satisfy the asymptotic identity 



Thus 
Explicitly, 



ri+S^#T 2 ~i+{^r)} 2 s^- (i4.8) 

L r=l M J r=l M 

Pi(«— g Wr)ftpi(r), ^(f)=| (A(r)} 2 -i wmW)- (14.9) 

^-15^-6 ?«. JL (14.10) 

48(* 2 -l) 3 96(* 2 -l)*f4 64f 3 



For the zeros of the derivatives of the parabolic cylinder functions, we may show by 
similar analysis that 



where 



4..-2» M r*(»+^+a^ + ..l (i4.il) 

L MM J 



3=M~%«, (14.12) 



&(»=*(», ft(f)=ft(r)*'(f), &(r)=A(«* / (f)+|{A(f)}*'*(f)»---. ( 14 - 13 ) 

and the &(£) are given by eq (7.22) of [8]. 
The equation corresponding to (14.5) is 

C/K< s )=±(2x)i|r(i-a)| i (-a-i<%)" i (- < ^)~ i , (14.14) 
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and from it we derive 

U(a,u' a .,)~ ( -^ {rQ+i/)}*Ai(a;)<»(/3) [l+g^]' (14-15) 

in which the coefficients satisfy the identity 

For the function U(a,z), the expansions corresponding to (14.1), (14.6), (14.11), and 
(14.15) are obtained on replacing the symbols U, u, Ai, and a s by U, u, Bi, and b 8 , respectively. 

(ii) a negative; other zeros. Since the zeros of the function Ai are all real and negative, it 
follows from (8.11) that for all sufficiently large positive values of the parameter /*, the zeros of 
C7(— |/x 2 , ju^V2) in the ^-domain T(0) (the unshaded region of fig. 7(a)), all lie in the interval 
— 1<^<0- This is also true of the zeros of the derived function U f {—^ji 2 ,fd^2). The func- 
tion U{— J/x 2 ,^V2) and its derivative however, have complex zeros in this region corresponding 
to the complex zeros of Bi and Bi' (cf. [8], appendix). 

In the complementary region CT(0) the distribution of the zeros of U(— §^ 2 ,^V2) can 
be investigated by reversion of the expansion (9.7). The situation is more complicated; the 
pattern of the zeros depends largely upon the non-integer part of \\£. It should be noted that 
expansions of zeros in the interval — \<Ct<C\ of fixed enumeration, counted away from t= — 1, 
come within this category. 

(hi) a complex. The extension of the results of this section to general values of a presents 
no difficulty. We readily see, for example, that when e <arg /x <ir— e the zeros of U(— §/z 2 , fd-^2) in 
the ^-domain T(arg/x) are asymptotically close to the principal curve AQ, illustrated in figure 4. 

15. Zeros of W(a,x) and W (a,x) 

We denote the sth positive zeros of the functions W(a,x), W'(a,x), W(a, — x) and W(a,— x), 
by Wa, s , v>a,s, w a ,s and w' at8 , respectively. 
(i) a positive. Writing n=j(2a) and 

a=n-%b s , /3=/T*W, (15.1) 

where b s and b f s are the sth negative zeros of Bi and Bi' respectively, we find by reversion of 

(12.18) and (12.20) 

L M M J 

2h^ a W'(a,w atS ) 2"M M *Bi' (&,)*(-«) [i+^(-y^z5fl (15.3) 

«\.~2*„ U(-fi-Jktfi+*tfi_ . . 1 (15 . 4) 

2le*-Pr(a,<.)~2M M -*Bi (b'M-fi) \l+±; (-)' ^=^1 (15.5) 

L r=l /* _J 

as a->+°°, uniformly with respect to s. Here ^(f), 0(f), p T (£), P r (f), <2>(D an d # r (f) are the 
functions defined in section 14 (i). 

With a, j8 given by (14.2) and (14.12) instead of (15.1) and Bi, b replaced by Ai, a, the 
right-hand sides of (15.2), (15.3), (15.4), and (15.5) are the corresponding expansions of w atS , 
—2-h~i ira W'(a,—w a , s ), Wa tS and 2~h~^ a W(a,— wi s ), respectively. 
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(ii) a negative. Let /x=V( — 2a), and 

*=r(I) (15.6) 

be the function inverse to f(£), defined by (11.7). Then by reversion of the expansion (12.28), 
we find 



W a 



08 ..~2*,[x+^+^+. • •] (15.7) 

as a->— oo , uniformly with respect to s. Here 

(15.8) 



•{^-}. 



and the leading coefficients are given by the formulas 

, , r 3 +6r , . 56r 9 +247r 7 +291r 5 -2520T 3 +3780r „ _ .. 

e 'W=24(r* + l)*' C2(r)= 5760(t*+1) s (15 - 9) 

The expansion of the associated value of the derivative can be obtained by substitution 

of (15.7) in the equation 

W'(a } w a , ,) = ±n4(^iV (15.10) 

obtained from eq (2.15) of [6] and (2.19). Thence we derive 

W f (a ) w a , 8 )^(-y2-i^(^+l)i[l+^^^] > (15.11) 

L r=l M J 



where 



„, , r 4 +15r 2 -6 - , , 112r 10 +489r 8 -98r 6 -34995r 4 +80460r 2 -7020 ,, B 1oN 

^ T ) = 48(r 2 + 1)3' ^ (T)=== " 23040(r 2 +l)6 " (15 ' 12) 



and 



In a similar manner we obtain from (12.29) 

J T i fi(r) , / 2 (t) , /icio\ 

'a.*~2 /x H s — | § — h . . . b (15.13) 

i mm i 



where now 



(15.14) 



, ,, r 3 -6r -,, 56r 9 -f247r 7 +651r 5 +2640r 3 -3420r /« * , e v 

/l(r)= 24(7+IF /■« — 5760(r 2 +l) 5 (15 ' 15) 



Finally, from 

^(a,<.) = ±T*Q<;-a)"*(^)"*, (15.16) 

we find 

W r (a,< s )~(-) s 2V*(r 2 +l)-iri+Z;^#T (15.17) 

where r is again given by (15.14), and 

„,, -t 4 -9t 2 +6 „,, -112t 10 -479t 8 -142t 6 +25725t 4 -67500t 2 +7380 ,._.„, 

^ l(r)= 48(r 2 +l) 3 ' ^M— 23040(r 2 +l)« " (15 " 18) 

If in the right-hand sides of (15.7) and (15.11) we take the value (15.14) for r instead 
of (15.8), we obtain the asymptotic expansions of w a>s and — W'(a, — w atS ), respectively. Sim- 
ilarly the right-hand sides of (15.13) and (15.17) with r given by (15.8) are the expansions 
of w'a tS and —W(a,—w' atS ). 
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It may be noticed that some of the numerical coefficients occurring in the expressions 
(14.4) and (14.10) for the coefficients in the reverted Airy- type expansions are the same as 
those in the expressions (15.9) and (15.12) pertaining to the reverted exponential-type expan- 
sions. This is no mere coincidence ; we may expect that by suitable re-expansion of the reverted 
Airy forms, and subsequent comparison with exponential-type expansions obtained by reversion 
of (5.11), we shall obtain explicit expressions for the coefficients p r and P r in terms of the 
e r and E T respectively, analogous to the formulas (8.12). Similarly for g r and Q r . 

16. Numerical Examples 

Example 1. To compute the smallest positive zero of W(l0,x) and the corresponding 
PF(10,a;). 

We use (15.2) and (15.3) with m=V20 = 4.47213 595, s=l. Entering table V of [3], 
we find 

6 a = — 1.17371 322, Bi'^) =0.60195 789. 
Hence 

f=-a=-AT^i = 0.15929 7392. 
The value of t is given by 

t(t 2 -l)^-\n{t+(t 2 -l)^}=~ f*=0.08477 1864. 

Solving this equation by successive approximation, using Newton's rule, we find 

*=1. 12489 6229=po(-o0, 
and thence from (14.4) 

^(-00 = 0.02960 04. 
Hence 

Wio,i~22X4.47213 595(1.12489 6229-0.00007 4001+ . . . ) = 7.11400 040. 

For the derivative, we compute from (14.7) and (14.10) 

lAK-a) = <K-a<)=0.88019 8053, P^ — a)=0.00942 2. 

Substitution in (15.3) then yields 

W'(lO,w iOA ) 2 ~ iT ^x0^9ml ^ t 1 - - 00002 356 + • • • ) = -1-39416 865X10-'. 

By inverse and direct interpolation in table III of [4], we find 

w 1M =7.11400 000, W , (10,i% u ) = -1.39416 86X10- 7 . 

The discrepancy in the values of w Wfl is undoubtedly due to neglect of the third term of the 
asymptotic series. 

Example 2. To compute the fifth positive zero of W f (—I0,x) and the corresponding 
W(-W,x). 

We use (15.13) and (15.17). From (15.14), (15.6), and (11.7) we see that the appropri- 
ate value of r satisfies 

r(r 2 +l)*+ln{r+(T 2 +l)i}=2f=2(5-i)7r/20 = 1.49225 6510. 
Solving by successive approximation, we find 

t=0.69390 0720, (t*+1)* = 1.21716 8110, 
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and thence from (15.15) 

/i(t) = -0.07269 508, / 2 (t) = -0.03321. 

Thus 

w / _ 10f6 ~2 i X4.47213 595(0.69390 0720-0.00018 1738-0.00000 0208+ . . . ) = 4.38746 276. 

From (15.18) we compute 

F t (r) =0.00919 198, F 2 (r)=^ 0.07882, 
and thence 

W(— 10y_ 10(5 ) 2 i (20)~ i (1.21716 8110)-*(1+0.00002 2980-0.00000 0493+ . . .) 

= - 0.50972 3297. 

The results of this example agree with values obtained by interpolation in table III of [4]. 
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Diffraction of electromagnetic waves by smooth 
obstacles for grazing angles, J. R. Wait and A. M. 
Conda, J. Research NBS 63D, No. 2, 18 (1959). 

The diffraction of electromagnetic waves by a convex cylin- 
drical surface is considered. Attention is confined primarily 
to the region near the light-shadow boundary. The complex- 
integral representation for the field is utilized to obtain a 
correction to the Kirchoff theory. Numerical results are 
presented which illustrate the influence of surface curvature 
and polarization on the diffraction pattern. Good agreement 
with the experimental results of Bachynski and Neugebauer 
is obtained. The effect of finite conductivity is also con- 
sidered. 

Path antenna gain in an exponential atmosphere, 

W. J. Hartman and R. E. Wilkerson, J. Research 
NBS 63D, No. 3, 273 {1959). 

The problem of determining path antenna gain is treated here 
in greater detail than previously. The method used here takes 
into account for the first time the exponential decrease of the 
gradient of refractive index with height, and a scattering cross 
section inversely proportional to the fifth power of the scatter- 
ing angle. Results are given for all combinations of beam- 
widths and path geometry, assuming that symmetrical beams 
are used on both ends of the path and that atmospheric 
turbulence is isotropic. The result appears as a function of 
both of the beamwidths, in addition to other parameters, and 
thus the loss in gain cannot be determined independently for 
the transmitting and receiving antennas. The values of the 
loss in gain are generally lower than the previous estimates for 
which a comparison is possible. 

Pattern synthesis for slotted-cylinder antennas, 

J. R. Wait and James Householder, J. Research NBS 
63D, No. 3,303 (1959). 

The radiation from a cylinder excited by an array of axial slots 
is discussed. A procedure for synthesizing a given radiation 
pattern is developed with particular attention being paid to a 
Tchebyscheff type pattern. Specifying the side lobe level and 
the width of the main beam, the required source distributions 
are computed for a number of cases. The effect of using a 
finite number of slot elements to approximate the continuous 
source distribution is also considered. 

Convexity of the field of a linear transformation, 

A. J. Goldman and M. Marcus, Can. Math. Bui. 2, 

No. 1, 15 (1959). 

The field F(A) of a linear transformation A of unitary n-space 
U n is the set of complex numbers z=(Ax,x), where x ranges 
over the unit sphere in U n . The convexity of F(A) (as a 
subset of the z-plane) is proved by a simple inductive argument 
which reduces the essential computations to the case n = 2. 

Graphical diagnosis of interlaboratory test results, 

W. J. Youden, Ind. Qual. Control XV, No. 11, 1 
(1959). 

This analysis of interlaboratory test results depends upon the 
availability from a number of laboratories of a single test 
result for each of two materials. The two results from each 
laboratory are used to plot a point using the a>axis for one 
material and the y-axis for the other. The resulting configura- 
tion of the points from the several laboratories permits deduc- 
tions on the prevalence and extent of laboratory bias, the 
presence of sampling variation and the occurrence of blunders. 
In addition the graph provides an estimate of the precision of 
the test procedure results. 



The calculation of the field in a homogeneous 
conductor with a wavy interface, J. R. Wait, Proc. 

IRE 47, 1155 (1959). 

The field at any depth in a homogeneous conductor with a 
wavy interface is calculated. For purposes of illustration, the 
field above the interface is taken to be a uniform plane wave 
traveling in the horizontal direction. 

Lower bounds for eigenvalues with application to 
the helium atom, N. W. Bazley, Proc. Nat. Acad. 
Sci. U.S. 45, No. 6, 850 (1959). " 

Let A be a self-adjoint operator with domain D in a Hilbert 
space. Suppose A = A fJ rA where A is self-adjoint and A' is 
positive definite. The eigenvalue problem for A, whose solu- 
tion we assume known, gives rough lower bounds. If u% 
(i= 1, . . ., k) are k discrete eigenvectors of A and if p t — (a')~ l 
Ui(i=l, . . ., k) exist then one can substantially improve the 
lower bounds. The theory is applied to helium atom operator. 
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